Astronomy & Astrophysics manuscript no. ohpaper'final © ESO 2013 

March 1,2013 



(N 

o 

(N 



(N 

<: 
O 

6 



oo 
a^ 

en 

o 



Letter to the Editor 

SOFIA observations of far-infrared hydroxyl emission toward 
classical ultracompact HII/OH maser regions 

T. Csengeri^ K. M. Menten', F. Wyrowski\ M. A. Requena-Torres\ R. Giisten', H. Wiesemeyer^ H.-W. Hiibers^'^, R 

Hartogh^, and K. Jacobs^ 

' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany e-mail: [ctimea;kment en] Ompi fr.de 

- Deutsches Zentrum fiir Luft- und Raumfahrt, Institut fur Planetenforschung, Rutherfordstral5e 2, 12489 Berlin, Germany 

^ Institut fiir Optik und Atomare Physik, Technische Universitat Berlin, HardenbergstraBe 36, 10623 Berlin, Germany 

* Max-Planck-Institut fiir Sonnensystemforschung, Max-Planck-StraBe 2, 37191 Katlenburg-Lindau, Germany 

' I. Physikalisches Institut, Universitt zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 



Received ... ; accepted 

ABSTRACT 

Context. The hydroxyl radical (OH) is found in various environments within the interstellar medium (ISM) of the Milky Way and 
external galaxies, mostly either in diffuse interstellar clouds or in the warm, dense environments of newly formed low-mass and high- 
mass stars, i.e, in the dense shells of compact and ultracompact HII regions (UCHIIRs). Until today, most studies of interstellar OH 
involved the molecule's radio wavelength hyperfine structure (hfs) transitions. These lines are generally not in local thermodynamical 
equilibrium (LTE) and either masing or over-cooling complicates their interpretation. In the past, observations of transitions between 
different rotational levels of OH, which are at far-infrared wavelengths, have suffered from limited spectral and angular resolution. 
Since these lines have critical densities many orders of magnitude higher than the radio wavelength ground state hfs lines and are 
emitted from levels with more than 100 K above the ground state, when observed in emission, they probe very dense and warm 
material. 

Aims. We aim to probe the warm and dense molecular material surrounding the UCHIIR/OH maser sources W3(0H), GlO.62-0.39 
and NGC 7538 IRSl by studying the ^Tli/2, J = 3/2 - 1/2 rotational transition of OH in emission and, toward the last source also the 
molecule's "Hj/t, J = 5/2 - 3/2 ground-state transition in absorption. 

Methods. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe these OH lines, which are near 1.84 
THz (163 ;um) and 2.51 THz (119.3 yum), with high angular (~16"/U") and spectral resolution (better than 1 km s"'). 
Results. We clearly detect the OH lines, some of which are blended with each other. Employing non-LTE radiative transfer calcula- 
tions we predict line intensities using models of a low OH abundance envelope versus a compact, high-abundance source correspond- 
ing to the origin of the radio OH lines. From the observed velocities and line-widths we can place constraints on the origin of the 
emission, and with detailed modeling we show for instance that the OH emission of W3(0H) comes from the UCHIIR and not from 
the envelope of the nearby hot-core. 

Conclusions. The far-IR lines of the OH molecule provide important information on the density and temperature structure of 
UCHIIRs. Taking a low-abundance envelope component is not sufficient to reproduce the spectra for W3(0H) and GIO. 62-0.39, 
a compact, high OH column density source - corresponding to the OH radio emitting sources is definitively required. 

Key words. Stars: formation - ISM: HII regions - ISM: molecules - Masers - Submillimeter: ISM 

1 . Introduction of the observed peculiar profiles with multiple, narrow velocity 

.^ components of polarized emission with highly anomalous 

^ . The hydroxyl radical (OH) was the first interstellar molecule relative hfs component intensity ratios, with t he 1665 MHz 

H ■ detected at radio wavelengths. After various searches, F = 2-1 line usually be ing the strongest (IWeinreb et all 

■ ^■' IWeinrebetal.1 (Il963h found absorption in the two strongest Il965l iBarrett & Roge rsl ll966 ^. The maser nature was clinched 

(F = 2 - 2 and 2-1, near 1667 and 1665 MHz, respectively) by interferometry, which showed that the OH emission near 

of the hyperfine structure (hfs) transitions within the molecule's the prominent HII region W3 had brightness temperatures 

^n3/2,7 = 3/2 ground state toward the supernova remnant Cas exceeding 10^ K and arose from compact (< 0''05) spots 

A. Subsequent studies found absorption toward all four of the spread over a few arc sec (thousands of AU) and significandy 

"18 c m lines" first toward the Galactic center ('Gardner et alj offset (by 14' or » 8 pc) fro m the radio c ontinuum emission 

119641) and, later, many galactic lines of sight (e.g., Goss 1968^ maxirntim of the HII re gion dRoeers et all [1966. MoraneTaTI 

Soon, also emission near 1665 MHz was discovered toward ll967l lDavies et alJl l967"). Most interestindv. the maser position 

prominent HII reg ions, i.e., sites of recent star formation turned out to be consistent with that of a newly detected radio 

JWeaver etaTIl 19651) . This and the detection of weaker emission sour ce much smaller tha n W3: the archetypical example of 

at the frequencies of the F - 2-2 and 1 - 2 OH hfs lines (at what lMezger et all (Il967h . based on their 2 cm radio continuum 

1667 and 1720 MHz) allowed l^inreb et alJ £1965) to identify observations, termed "a new class of compact H II regions 

of OH as the carrier of this peculiar emission. It was readily associated with OH emission sources ". This source, W3(OH) 

identified as due to maser action by the latter authors because is characterized by its youth (2300 yr, iKawamura & MassonI 
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Il998h . compa ctness (~ 0.02 pc) an d high electron densities 
(~ 2 10^ cm'l lDreher & Welchlll98lh and is still surrounded by 
an ex panding dusty, molecular envelope (e.g., Wyrowski et al. 
Il999l) . It is the archetypical ultracompact HII region (UCHIIR) 
and is discussed in this Letter along with the UCHIIRs 
GlO.62-0.39 and NGC 7538 IRSl. 

In addition to the lines from the ground state, transitions be- 
tween hfs levels of various rotationally excited states with ro- 
tational quantum number J from 5/2 to 9/2 from the lower 
energy ^113/2 ladder and from 1/2 to 5/2 in the ^ni/2 ladder 
hav e been detected at radio frequ encies between 4 and 23 GHz 
(see lCesaroni & Walmslevl[l99lL for a summary of the observa- 
tions and modeling). These lines probe levels with energies up 
to 5 1 1 K above ground, whose populations are determined by an 
interplay of strong far-infrared (FIR) radiation from warm dust 
and collisions in dense gas since the FIR rotational lines that 
(de)populate them, at frequencies > 1.8 THz, have high critical 
densities. Radiative transfer calculations show that many of them 
are not in local thermodynamic equilibrium (LTE), which gives 
rise to maser emission in some lines and enhanced absorption, 
caused by "overcooling", in others, see Cesaroni & Walmslev 
( Il991h . As these authors show, just the fact that some lines do 
show maser emission and others do not places interesting con- 
straints on kinetic temperature and density. 

The modeling of the various cm OH hfs lines delivers a 
consistent picture of the OH in warm, dense gas. Nevertheless, 
for a full characterization, observations of the rotational OH 
transitions connecting hfs levels from different rotationally ex- 
cited states are desirable. The frequencies of these transitions, 
within the ^113/2 and 1/2 ladders, and also between them, 
all lie in the supra -THz range. Using the Kuiper Airborne 



I 



Observatory (KAO), IStorev et all (Il981h made the first detec- 
tion of the ^113/2,7 = 5/2 - 3/2 ground-state transitions near 
2.51 THz (119. 3 yum) in absorption toward Sgr B2 and in emis- 
sion from the shocked gas in the Kleinmann-Low (KL) neb- 
ula in the Orion molecular could 1. Detections of other OH 
lines followed, including those of the other lines discussed 
in the present Letter, the ^ni/2,7 = 3 /2 - 1/2 transitions 



near 1.83 THz (164 yum), also in Orion (iViscuso et al. 19851). 
The Infrared Observatory (ISO) brought a significant improve- 
ment in sensitivity, although neither in angular nor in spec- 
tral resolution. Again, studies concentrated on the exception- 
ally bright and extended sources Sgr B2 and Orion-KL. Multi- 
transition data aff'orded by ISO's long wavelength spectrome- 
ter allowed detailed excitation studie s of both sources in several 
hydroxyl isotopologes ('^'^^ ' ^OH) ( jGoicoechea & Ce rnicharol 
[2002; Goicoechea et al] |2006l: [Poleha mpton et al. 2003). 

Herschel and SOFIA have brought vast improvements over 
both the KAO and ISO in terms of sensitivity and angular 
and velocity resolution. Both are flying heterodyne spectrom- 
eters affording bandwidths of several hundred km s ' even at 
THz frequencies and velocity resolutions better than 1 km s ', 
namely the Heterodyn e Instrument for the Far Infrared (HIFI, 
Ide Graauw et al.ll2010h and the Germ an Receiver for Astro nomy 
at Terahertz frequencies (GREAlQ, H evminck et alJl2012h . The 
increased sensitivity allows studies of more compact regions, 
namely the above mentioned dense molecular envelopes asso- 
ciated with UCHIIRs and OH masers. Recentlv. lWampfler et alJ 
(l2011h presented HIFI observations of the ^111/2, 7 = 3/2-1/2 



' German Receiver for Astronomy at Terahertz frequencies. 
GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 



triplet of OH hfs transitions at 1837.8 GHz (163.1 ^im) toward 
the high-mass star-forming region W3 IRS 5, for the first time 
resolving the hfs structure. A multi-rotational OH line study 
of the Orion Bar photon-dominated region wi th the R^.CS in- 
strume nt onboard Herschel was presented by .Goicoechea et al.l 
(l2011h . Here we report in Sec. |2] our SOFIA observations of 
emission in this and the other A-doublet component at 1834.7 
GHz (163.4 jum) toward the well-known UCHIIRs/OH maser 
sources W3(OH), GlO.62-0.39 and NGC7538. The latter two 
have similar warm, dense molecula r environmen ts, and also 
sho w radio OH mase r emission (Argo n et al.ll2000l) and absorp- 
tion (Walmslev et al. 1986). Toward the latter source we also ob- 
served absorption from the ^113/2,7 = 5/2 - 3/2 ground-state 
transitions near 2.51 THz (1 19.3 jum). 

2. Observation and data reduction 

The OH ^ni/2, J - 3/2-1/2 rotational lines were observed with 
the Pl-instrument GREAT onboard SOFIA. The observations 
were caiiied out during flights no. SS02-OCF4/03 and BS02- 
01/06/07. The signal was detected at 1837.8 GHz (163.1 yum) 
using the L2 band of GREAT, while the pair of transitions from 
the other A-doublet at 1834.7 GHz (163.4 yum) lies in the image 
sideband. The beam is 15:'5 at 1837.8 GHz. The AFFT back- 
end was used with 1.5 GHz bandwidth and 212 kHz spectral 
resolution. The average system temperature is 4700-5300 K. 
Toward NGC 7538 IRSl the ^n3/2,J = 5/2" - 3/2+ rotational 
line at 2514 GHz was observed using the M-band receiver See 
Wiesemeyer et al. (2012) for details of these observations. The 
pointing was established with the optical guide cameras (with an 

accuracy of 5")- 

The data were calibrated using standard procedur es based 
on the KOSMA/GREAT calibrator (iGuan et al.ll2012l) . for the 
data analysis the GILDASQ software was used. The forward effi- 
ciency was set to 95% and the data were converted to a Rayleigh- 
Jeans equivalent T^b.Rj scale using a beam efficiency of 5 1 % for 
the L2 and 58% for the M band, respectively. 

The spectra were first summed up to define windows ex- 
cluding the line signal for the baseline fitting. Then a third or- 
der polynomial baseline was fitted and subtracted from each 
scan individually. The ^ni/2,7 - 3/2 - 1/2 spectra were av- 
eraged with noise weighting and are smoothed to ~1.2 km s"^ . 
We reach on average ~0. 3-0.4 K noise level per channel. The 
^n3/2,7 = 5/2" - 3/2+ line was reduced in a similar fashion. 
Here, the spectrum is smoothed to 0.9 km s ' velocity resolu- 
tion and has a 0.4 K noise level. 



3. Results 

The A-doublet lines of the OH ^111/2, 7 = 3/2-1/2 rotational 
transitions are detected, and as is characteristic of this transition, 
they are seen in emission (Fig.[T]i toward all three high-mass star- 
forming sites. The F = 2+ - 1" and F = 1+ - 0" hfs components 
(at 1837.8168 and 1837.8370 GHz, respectively) are separated 
by 3.2 km s ' and are therefore spectrally marginally resolved. 
The more separated and weaker F = 1+ - 1" component (at 
1837.7466 GHz) would be spectrally resolved, but lies within 
the noise level toward W3(OH) and NGC 7538 IRSl. Because 
GREAT is a double-sideband receiver, the other component 
of the A-doublet (a blended triplet at 1834.7355, 1834.7474, 
1834.7504 GHz, for the F= V - l+,2- - 1+, V - 0+ transi- 
tions, respectively) could be recorded in the image band. For 



See ht tp://www.iram.fr/IRAMFTVGILDAS7l 
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Table 1. Summary of the derived line parameters of the ^ni/2, J - 1^17. - 1/2 transition 



Source 



Position T„,,,j;j 

RA[J2000] Dec[J2000] [K] 



Vlsr 

km s" 



Av 
km s" 



total r 



[K] 



W3(0H) 

GlO.62-0.39 

NGC7538IRS1 



02:27:03.90 
18:10:28.64 
23:13:45.36 



61:52:24.6 
-19:55:49.5 
61:28:10.5 



1.83: 
1.34: 
1.04: 



:0.34 
:0.29 
:0.34 



-45.70 
-3.17 = 
-57.80 



t0.31 

:0.51 

tO.43 



7.54 ±0.87 
9.50 ±1.15 
5.46 ±1.00 



0.1 
0.1 
0.1 



2 40.2-5.1 
5 30.2 - 3.7 
5 24.1-3.5 



Notes. Line intensities are derived for the sum of the three hfs lines in each A-doublet of the 'ni/2, J = 3/2 - 1/2 transition. We give the range of 
total r values that give fits within the noise level and the corresponding T„ values. The total t corresponds to the sum of optical depths of all hfs 
components in both doublet components. 



GlO.62-0.39 the F = 1+ - 1" line partially overlaps with the 
emission from the image band. 

For all sources, the observed line profiles of this transition 
exclude the presence of a broad outflow component within our 
noise level of ~0.4 K. The lines are turbulence-dominated and 
well-fitted with Gaussians within this noise level. The derived 
line -widths are consist ent with those of other high-density trac- 
ers (iPlume et al.lll997h . suggesting an origin from a dense, tur- 
bulent medium. 

Toward NGC 7538 IRSl, the blended 2514 GHz lines of 
the ^113/2,7 = 5/2" - 3/2^ rotational transition were also ob- 
served and are clearly seen in absorption. The tiiplet of lines 
from the other component of the A-doublet at 2510 GHz was 
not observed. The level of continuum emission was unstable and 
could not be determined reliably. 

Given the OH molecules's substantial dipole moment 
(1.66 D), these high-frequency rotational transitions have high 
ciitical densities > 10^ cm"^. As a consequence, the LTE as- 
sumption is unlikely to apply even in the usual high-density en- 
vironment of massive protostars, and FIR radiative pumping is 
expected to play a role in the excitation. Therefore ideally one 
should fit the individual hfs components to derive their inten- 
sity ratio, however, the presented spectra do not have sufficient 
signal-to-noise ratio for that. We used, hence, the HFS fitting 
method with CLASS and simultaneously fitted all six hfs com- 
ponents of the doublets with fixed optical depth, t. Best fits, with 
a residual rms within the noise, are summarized in Table[T] The 
spectra's quality is moderate. For example, the broad feature un- 
derlying the 1837 GHz line toward W3(OH) is caused by an im- 
perfect spectral baseline. 

4. Analysis: Radiative transfer modeling 

We used RATRAN ('Hogerh eiide & van der I^l2000h to model 
OH emission from protostellar envelopes, where the effect of 
dust heating was included by adopting model s from the literature 
dvan der Tak et al.ll2000l: iMueller et al.ll2002 '). Only for W3(OH) 
we present the ^Il\/2,J - 3/2 - 1/2 line intensities calculated 
for a compact (few arcsecond size) source associated with OH 
maser emission with the larg e velocity gradient (LVG) model of 
ICesaroni & WalmslevI ( 1199 lb . In addition to the far-IR pumping, 
this takes into account line overlap effects and the hfs structure 
of the lines to give a detailed prediction for the cm radio lines of 
OH. 



4.1. W3(0H): OH emission dominated by fhe UCHII region 

ICesaroni & WalmslevI ( 1199 lb created an LVG codejj to model the 
line properties (emission, absorption and population inversion) 
for the cm radio lines of OH toward W3(OH), where the 4 GHz 
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Fig. 1. Upper three panels: the ^Yli/2,J - 3/2 - 1/2 rotational 
lines of the three UCHIIRs in emission. The gray dashed line 
shows Vlsr- Red lines indicate the position of the hfs components 
of the - parity line of the doublet, while blue labels show the hfs 
of the + transition. (Labels are shifted for better visibility.) The 
dashed green line shows the t - O.l hfs fit to the spectra. The 
lower sideband rest velocity is shifted according to the v;,, of 
the individual sources and is indicated on the upper axis of each 
plot, while the rest velocity scale is the same for all sources. The 
lowest panel shows the ^1X3/2,7 = 3/2 - 5/2 transition toward 
NGC 7538 IRSl. 



This code uses the collisional cross-sections from lFloweJ (1 1989h . 
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hfs lines are seen in both emission and absorption. They show 
that this particular combination originates from a very narrow 
range of physical conditions. As a first step we reproduced their 
model and looked at the predictions for the far-IR lines. 



The ICesaroni & WalmslevI (Il99lh model describes a small 
(few arcsecond) compact so urce based on interferometric obser- 
vation s from the literature (iGuilloteau et al.l[l985l IWilson et aP 
I1990I) . They qualitatively reproduce the radio line observations 
using Tgas - Tdusi = 151 K, including the presence of both inter- 
nal dust and line overlap effects. The model uses an abundance, 
Xou = 2 X 10"^ and a velocity gradient of 200 km s"' pc"' 
(see also ApplAli. Line overlap is taken into account and the ve- 
locity range available for non-local line overlap is fixed in this 
model to AV=0.7 km s ' . We reproduced their model and looked 
into the far-IR lines (Fig. lA.ll upper panel). The observations 
agree well with the predictions of the model in a density range 
of hhj = 2-3x10^ cm"^ with a corresponding r^v ~ 50-60 K. 
We note, however, that the model largely overpredicts the opti- 
cal depth. Using this LVG code we also tested the importance of 
internal heating by dust, i.e. far-IR radiative pumping from the 
dust continuum photons. As shown in Fig. lA.lK lower panel) we 
find that models without a radiation field strongly underestimate 
the line intensity by an order of magnitude. 

Because the kno wn hot-core source, W 3(H20), is only 6" 
offset from W3(OH) dWvrowski et al.ll 19991) . we tested whether 
adding an envelope component could influence our model re- 
sults. We used RATRAN (Hogerh eiide & van der Tak 2000) and 
adopted an envelope model from van der Tak et al.l (l2000h . with 



a power-law density profile n(r)=no(r/r())'' (no = 5.3x10 cm" , 
p = -2) and a temperature gradient with an exponent of -0.4, 
adding up to a total luminosity of 2x10"* Lq. We adopted an abun- 
dance of YoH = 0.8 X 10"*^, which is observed in a similar en- 
vironment by IWampfler et al.l (1201 ll) . These models predict very 
low line intensities (T/, ~ 0.3 K), therefore we conclude that 
for W3(OH) such an envelope model cannot significantly con- 
tribute to the observed line intensities. This result is consistent 
with the observed v/^^ ~ -46 km s"' of the emission, which 
coincide with the UCHII region, W3(OH), rather than the neigh- 
boring hot-core, W3(H20). 



4.2. NGC 7538 IRS1: OH emission dominated by an 
envelope component 

For this source, the 2514 GHz ^113/2,7 = 5/2 - 3/2 transition 
was also observed. The line is seen in absorption and although 
the absorption feature seems to be saturated, the continuum level 
is uncertain. Hence, the only constraint we can deduce is that 
it appears in absorption. (The similar line-width of the hfs fits 
suggests that the optical depth may be similar to that of the 
^111/2, 7 = 3/2-1/2 transition.) 

We used RATRAN to model an envelope with a power-law 
density and temperature profile with no=5.3xl O"^cm~^, p - -1-0 
and a t otal luminosity of 1.3xlO^L0 following , van der Tak et al.l 
(I2OOOI) . Assuming a constant abundance of ;ifoH = 0.8 X 10"^ 
we reproduced the observed emission line intensities and the ab- 
sorption feature as well (Fig. IB.lt . This indicates that an enve- 
lope model can explain rotational transitions of the emitting and 
absorbing OH gas. 



4.3. G10.62-0.39: OH emission from an envelope and a 
compact source 

Similarly as for the other sources, we produced models with 
RATRAN, where the physical model of the envelope is based on 
models of iPlume et alj d 19971) and lMuelleret all (|2002|) . adopt- 
ing a power-law density profile with an exponent p - -2.5, 
no=1.2xlO' cm"^ and a total luminosity of 9.2xlO^L0. These 
calculations produce a narrower emission feature than observed 
in our OH spectrum and underestimate the line intensity, giv- 
ing Tl ~ 0.6 K. Either a higher OH abundance or the addition 
of a compact component corresponding to the UCHII region is 
needed. 

5. Conclusions 

We have detected the OH ^ni/2, J - 3/2 - 1/2 rotational lines 
in emission toward three classical UCHIIRs and the ^113/2, J = 
3/2 - 5/2 rotational line in absorption toward one of them, 
NGC 7538 IRSl. Our modehng shows that the far-IR radia- 
tion field plays an important role in the excitation conditions. 
Except for NGC 7538 IRSl, an envelope model with a low OH 
abundance (yoh - 0.8 x 10"*^) cannot explain the observed line 
intensities alone; an additional compact component with high 
density and high OH abundance is needed. We have shown for 
W3(OH) that the emission is dominated by the UCHIIR and not 
the nearby hot-core. Given the particular excitation properties 
of the radio lines of the OH molecule, a detailed modeling of 
their properties (emission, absorption, or population inversion) 
and observations of more rotational transitions of OH can pro- 
vide more constraints on the physical conditions. 
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Appendix A: LVG model of W3(0H) 

We rep roduced the model of W3(OH) of ICesaroni & WalmslevI 
(Il99lh . which qualitatively models the radio lines of OH, and 
show here the predictions of this model for the 1834 and 1837 
GHz lines. We adopted their parameters, assuming a source size 
of 2" (0.01 pc), a velocity gradient of 200 km s"' pc"' and an 
OH column density between 2x10'"* cm"^ to 2x10"' cm"^. To 
compare these intensities with the observed values we used a 
beam size of 15''5 for SOFIA. In Fig. lA.ll two models are shown, 
with and without an internal radiation field introduced by the 
present of warm dust. As shown, the latter model underestimates 
the observed line intensities. 

This code is well-adapted for the OH molecule because it 
includes line-overlap effects, however, for a detailed modeling 
the radio OH lines need to be taken into account. To perform this 
modeling for all the sources is beyond the scope of this letter. 
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Fig.A .l. Results of the LVG models of ICesaroni & WalmslevI 
(Il99lh showing only the predictions of excitation temperature 
and hne intensity for the 1834 and 1837 GHz lines. The gray 
dashed area indicates the observed line intensity including the er- 
rors. Light green and blue lines correspond to the sum of the hfs 
lines in each component of the doublet. Upper panels: Models 
with internal radiation field. Lower panels: Models without in- 
ternal radiation field. 
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Fig.B.l. Our RATRAN model of the 1834 GHz and 2514 GHz 
OH line in emission and in absorption, respectively. Black line 
shows the observed spectra with baselines removed, the red line 
shows the model. 



Appendix B: RATRAN model of NGC7538 IRS 1 

RATRAN is based on a Mo nte-Carlo approach to perform ra - 
diative transfer calculations (iHoeerheiide & van der Takll2000h . 
It is therefore well-suited for detailed modeling of protostellar 
envelopes. A radiation field from dust heating was implemented 
in the models. Note that as shown in Fig. lB.ll a line profile can 
be produced. 



